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Abstract 

We construct an approximate solution which describes the gravitational emission from a naked 
singularity formed by the gravitational collapse of a cylindrical thick shell composed of dust. The 
assumed situation is that the collapsing speed of the dust is very large. In this situation, the metric 
variables are obtained approximately by a kind of linear perturbation analysis in the background 
Morgan solution which describes the motion of cylindrical null dust. The most important problem 
in this study is what boundary conditions for metric and matter variables should be imposed at 
the naked singularity. We find a boundary condition that all the metric and matter variables are 
everywhere finite at least up to the first order approximation. This implies that the spacetime 
singularity formed by this high-speed dust collapse is very similar to that formed by the null dust 
and thus the gravitational emission from a naked singularity formed by the cylindrical dust collapse 
can be gentle. 
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I. INTRODUCTION 

The cosmic censorship hypothesis is a crucial ansatz for the theorems of black holes[J|. 
There are two versions for this hypothesis. For spacetimes with physically reasonable matter 
fields, the weak version claims that a spacetime singularity as a result of generic non-singular 
initial data is not visible from infinity, while the strong version claims that a spacetime sin- 
gularity developed from non-singular initial data is invisible for any observer. A singularity 
censored by the strong version is called a naked singularity, while a singularity censored by 
the weak version is called a globally naked singularity. However, previous theoretical stud- 
ies have revealed several candidates for the counter-examples of this hypothesis^, U, [3, 
Although further detailed and careful studies about these candidates are necessary, these 
might have some physical importance. 

If a globally naked singularity forms, what can we detect from there? In connection with 
this issue, Nakamura, Shibata and one of the present authors (KN) have proposed a conjec- 
ture; if globally naked singularities form, large spacetime curvatures in the neighborhood of 
these can propagate away to infinity in the form of gravitational radiation, and as a result, 
almost all of the mass of these naked singularities is lost through this large gravitational 
emission J]]]- If this conjecture is true, the formation processes of the globally naked singu- 
larities might be very important for gravitational- wave astronomy ^, Q| . The cylindrically 
symmetric system will play an important role in getting significant information about this 
issue by the following reasons; in the asymptotically flat case, the singularities are neces- 



sarily naked jloL lllj: there is a degree of gravitational radiation; this system is very simple. 



There are a few numerical studies on gravitational emissions by cylindrically symmetric 



gravitational collapse 
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Recently, the present authors have investigated the high-speed collapse of cylindri cally 
symmetric thick shell composed of dust [15] and perfect fluid with non- vanishing pressure jl6|. 
In these studies, we use the high-speed approximation which is a kind of the perturbation 
analysis in the background Morgan solution describing the motion of a cylindrical null dust. 
In the first paper, it has been revealed that in the gravitational collapse of a cylindrically 
symmetric thick dust shell, the thinner shell leads to the larger amount of gravitational 
radiation. This result gives a resolution of an apparent inconsistency between the previous 
works. In the second paper, we showed that the pressure decelerates the collapsing velocity 
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so significantly that the high-speed approximation scheme breaks down before a singularity 
forms if the equation of state is moderately hard. For example, in the case of mono-atomic 
ideal gas which seems to be physically reasonable in the high-energy state due to the asymp- 
totic freedom of elementary interactions, the high-speed collapse is prevented by the pressure 
and thus there is a possibility that the pressure halts the formation of a spacetime singularity 
in the case of this cylindrical gravitational collapse. However, we should note that if initial 
collapsing velocity is very large, a region with very large mass concentration can be realized 
in the neighborhood of the symmetric axis r = 0. In this region, the tidal force for free 
falling observers can be so large that general relativity breaks down and the quantum the- 
ory of gravity will be necessary to understand the physical processes realized there. Such a 



region can be regarded as a singularity for general relativity, called a 'spacetime border' [17 1. 
By contrast, in case that the equation of state is sufficiently soft, the high-speed collapse is 
maintained until a globally naked singularity forms. 

In this paper, we investigate the gravitational waves from a naked singularity formed 
by the high-speed gravitational collapse of a cylindrical thick shell composed of dust. In 
the previous paper pjjj, we have investigated the generation process of those in the causal 
past of the Cauchy horizon associated with the naked singularity, i.e., the region which does 
not suffer the influence of the naked singularity (for details of the Cauchy horizon, causal 
future and causal past, see, for example, the textbook by Wald[18j). On the other hand, 
in this paper, we focus on the generation process of the gravitational waves in the causal 
future of the naked singularity. In order to know what happens in the causal future of the 
naked singularity, we need specify the boundary condition at the naked singularity. This is 
equivalent to defining the naked singularity as a physical entity. We would like to show that 
in the high-speed approximation scheme, there is a boundary condition that all the metric 
and matter variables are everywhere finite at least up to the first order. There are several 



studies on how to fix the boundary condition 



in static naked singular spacetimesj 



B H 21, 22, Q 



or test fields or gravitational perturbations 



By contrast, the present case is 



dynamical, and as far as we are aware, this is the first example of gravitational emissions 
from a naked singularity developed from non-singular initial data. 

This paper is organized as follows. In Sec. II, we present the basic equations of the cylin- 
drically symmetric dust system. In Sec. Ill, we derive the basic equations in the high-speed 
approximation scheme which describe the gravitational collapse with velocity almost equal 
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to the speed of light c. In Sec. IV, we consider the boundary condition at the background 
naked singularity and then present solutions for the basic equations given in Sec. III. In 
Sec.V, we study the gravitational emissions from the naked singularity. Finally Sec. VI is 
devoted to summary and discussion. 

In this paper, we adopt c = 1 unit and basically follow the convention of the Riemann 
and metric tensors and the notation in the textbook by Wald 18]. 

II. CYLINDRICALLY SYMMETRIC DUST SYSTEM 

In this paper, we focus on the spacetime with whole-cylinder symmetry which is defined 



by the following line element 



24] 



ds 2 = e 2 ^ (-dt 2 + dr 2 ) + e- 2 ^R 2 dif 2 + e 2 ^dz 2 . (2.1) 
Then Einstein equations are 

7 ' = - i? 2 )" 1 S^RR' (ip 2 + ^ - 2RRipijj' + R'R" - RR' 

- 8nGV=g (R'T\ + RT r t )}, (2.2) 
7 = - (r' 2 - i? 2 ) _1 jiLR (V> 2 + ?// 2 ) - 2RR'i)i)' + RR" - R'R' 

- 8ttG^9 (RT\ + R'T r t )\, (2.3) 

7-7" = /-^ 2 -^TV (2.4) 
R-R" = -8ttGv^ (T\ + T r r ) , (2.5) 

$ + f ^ - r ~ f V = -^v^ (T\ + T\ - T z z + T%) , (2.6) 

where a dot means the derivative with respect to t while a prime means the derivative with 
respect to r. 

As mentioned, we consider the dust fluid as a matter field. The stress-energy tensor is 

T^u = pu^Uu, (2.7) 

where p is rest mass density and u M is a 4- velocity of the fluid element. Due to the assumption 
of the whole-cylinder symmetry, the components of the 4- velocity are written as 

= u t (l,-l + V, 0, 0). (2.8) 
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By the normalization u^u^ = — 1, u l is expressed as 



«* = . = (2.9) 
'V (2 - V) 



Here denoting the determinant of the metric tensor g^ v by g, we introduce new density 
variable D defined by 

D:= = R f~*'\. (2.10) 

y]V(2- V) V(2-V) 

The components of the stress-energy tensor are then expressed as 

^^T\ = -e^D, (2.11) 
y/=g~T r t = e^'Dil -V) = -^gT\, (2.12) 
yf^T% = e 7 "^ (1 - V) 2 D, (2.13) 

and the other components vanish. 

The equation of motion V a T a ^ = leads to 

d u D = ~(DV)' + | (1 - V) {2d u (4> - 7 ) - V(j> - 7)} , (2.14) 

Dd u V = (1 - y)0 tt Z> + ~ {V(l - V)£>}' - ^ {29 u (^ - 7) - V{iP - 7)} , (2.15) 

where u = t — r is the retarded time and d u is the partial derivative of u with fixed advanced 
time v = t + r. The first equation comes from the t-component while the second one comes 
from the r-component. The z- and ^-components are trivial. 

To estimate the energy flux of the grav itational radiation, we investigate C-energy E 
and its flux vector proposed by Thorne[24j. C-energy E = E(t,r) is the energy per unit 
coordinate length along z-direction within the radius r at time t, which is defined by 

E = ^{l + e- 2 ^(R 2 -R' 2 )}. (2.16) 

The energy flux vector J M associated to the C-energy is defined by 

^ff'-i' ' )' (217) 

By its definition, J M is divergence free. Using the equations of motion for the metric variables, 
we obtain the following expression of C-energy flux vector, 

~gJ l = ^{ RR '^ 2 + ^) - 2RR>H/ - SnG^iR'T't + RT r t )Y (2.18) 

—gf = t^l^RR^ 2 + <// 2 ) _ 2RR'iPiP' - 8nG^(R'T r t - i?T r r )j, (2.19) 



and the other components vanish. 



III. HIGH SPEED APPROXIMATION SCHEME 

Let us consider the ingoing null limit of the cylindrical dust fluid. From Eqs. (j2.11j) - (j2.13|) . 

the stress-energy tensor is written in the following form, 

e 3W-7)£) 
~R~ 



T^u — — k^k u , (3-1) 



where 

fc" = (1,-1 + V, 0, 0). (3.2) 

The timelike vector k^ becomes the ingoing null vector in the limit of V — > 0+. Hence in 
the limit of V — > + with D fixed, the stress-energy tensor agrees with that of the collapsing 
null dust. This means that in the case of very large collapsing velocity, i.e., < V <C 1, 
the dust fluid system will be well approximated by the null dust system. Then we treat the 
"deviation V of the 4- velocity from null" as a perturbation and perform linear perturbation 
analyses. 

In the case of the complete null dust V = 0, the solution is easily obtained as 

-4> = 0, (3.3) 

7 = 7b0), (3.4) 

R = r, (3.5) 

87rG.De 7 = (3.6) 

where 7r(u ) is an arbitrary function of the advanced time v. This solution was first obtained 
by Morgan 2^ and was studied subsequently by Letelier and Wang and NolanJ^ in 
detail. Assuming the null energy condition, the density variable D is non-negative. This 
means that 7b is non-decreasing function of v. 

The situation of the collapsing null dust solution can be understood by Fig.l. In this 
paper, we assume that the density variable D has a compact support < v < v w which 
is depicted by a shaded region in Fig.l. We can find in Eqs. (|3.1|) and (|3.5|) that if D does 
not vanish at the symmetric axis r = 0, components of the stress-energy tensor T^ v with 
respect to the coordinate basis diverge there, and the same is true for the Ricci tensor by 
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t 




FIG. 1: Morgan's cylindrical null dust solution. There is the null dust in the shaded region; 
D(v) > for < v < v w . The dashed line on t-axis corresponds to the intermediate singularity 
at which an observer suffers infinite tidal force although any scalar polynomials of the Riemann 
tensor do not vanish. On the other hand, the dotted line on t-axis is the conical singularity. The 
Cauchy horizon is short dashed line of t = r. We estimate C-energy transfered through the null 
hypersurface v = v w from u = to u = v w , which is depicted by dot-dashed line. 

Einstein equations. This is a naked singularity which is depicted by the dashed line at r = 
in < t < v w in Fig.l. Although all the scalar polynomials of Riemann tensor vanish 
there, freely falling observers suffer infinite tidal force atr = 0in0<t< v w ; this type of 
singularity is called p.p. curvature sing 

ularityH- 

The Cauchy horizon associated with this 
intermediate naked singularity is represented by the short dashed line t = r in Fig.l. On 
the other hand, the region at r = in t > v w is conical singularity which is depicted by the 
thick dotted line in this figure. 

We introduce a small parameter e and assume the order of the variables as V = 0(e) and 
if) = 0(e). Further we rewrite the variables 7, R and D as 

e 7 = e ^(l + 5 7 ), (3.7) 
R = r(l + 5 R ), (3.8) 
D = D B (1 + 6 D ), (3.9) 
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and assume that 5 7 , 5 R and 5d are 0(e), where 

D " = sd^iB- < 310 > 

We call the perturbative analysis with respect to this small parameter e the high speed 
approximation scheme. 

The 1st order equations with respect to e are given as follows; the Einstein equations 
(Q-tESD lead to 

<V = 8nGD B e^ {5 7 - ^ + - 2d v (r5 R )} + (r8 R )", (3.11) 
5 7 = 8nGD B e^ {5 7 - i/j + 5 D - 2d v (r5 R ) - V} + (r5 R )', (3.12) 
<5 7 -5 7 " = 0, (3.13) 
r5 R - (r5 R )" = lQ7iGe^ B D B V, (3.14) 
# - - V = — c^DbV; (3.15) 

the conservation law 1)2.14)1 leads to 

^ D + S ^^-±-^V-\(v-^v); (3,6) 
the Euler equation ()2.15)) becomes as 

d u V = 0, (3.17) 

where we have used Eq. (j3.16|) . 

C-energy E up to the first order of e is given by 

1 



E = - 



1 - e- 27B + 2e" 27B {<5 7 - (rS R )'}] . (3.18) 



From Eq. (j3.6|) . we can easily see that 7b is constant in the vacuum region D B = 0. Further, 
from Eqs. ()3.11j) and ()3.12|) . we find that <5 7 — (r5 R )' is also constant in the vacuum region. 
Therefore up to the first order, C-energy is constant in the vacuum region. This means that 
in the vacuum region, C-energy flux vector J M vanishes up to the first order and thus it is 
a second order quantity which is given by 

r 



C-energy flux vector takes very similar form to that of massless Klein-Gordon field. 
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IV. SOLUTIONS FOR DUST COLLAPSE 



In this section, we study the behaviors of perturbation variables in the causal future of 
the background naked singularity. Einstein equations are hyperbolic differential equations 
which send out the information of the background naked singularity into its causal future. 
Therefore in order to get the solutions of the basic equations for the perturbation variables 
in the causal future of the background naked singularity, we need specify the boundary 
condition for these at the background naked singularity. 

In general, it is very difficult to control the solutions by the boundary conditions due to 
the nonlinearity of Einstein equations. However, fortunately, there is a great advantage in 
the present situation; the perturbation variables can be regarded as test fields with sources in 
the fixed background Morgan spacetime. This seems to be similar situation to the problems 
investigated by Wald and a several researchers, i.e., the behavior of test fields in naked 
8mgular static spaeetitnesfl H H Q Q. However, there is a secant difference 
between the previous studies and the present case. In the previous cases, the effect of the 
naked singularity appears in the differential operator in the equations of motion for the test 
fields. On the other hand, in the present case, the manifestation of the naked singularity 
appears as source terms in the differential equations for the perturbation variables. Therefore 
we have to determine the boundary condition from different point of view from previous 
works. 

We do not know the law of physics at the spacetime singularity, i.e., the quantum theory 
of gravity and hence, at present, it is difficult to impose a physically satisfactory boundary 
condition at the naked singularity. There is a possibility that the quantum gravity puts 
the additional source terms on the right hand side of the differential equations (j3.11|) - ()3.15|) 
at the background naked singularity; the additional source term might be a delta-function 
like one causing gravitational wave bursts, or by contrast, might completely cancel out 
the original source term and, as a result, lead no gravitational radiation from the naked 
singularity. However, at first glance, it is expected that the geometrical structure of the 
spacetime investigated here is not so different from the background Morgan spacetime even 
in the neighborhood of the background naked singularity. Thus in order to search for such 
a solution, we require that the perturbation variables are expressed in the form of the 
Maclaurin series around the naked singularity. As will be shown below, this condition 
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guarantees that all the perturbation variables are everywhere finite up to the first order, 
and this is just what we expect. The cylindrical thick shell composed of dust also forms a 
naked singularity at the same coordinate position as the background naked singularity The 
Cauchy horizon associated with this naked singularity is located also at the same coordinate 
position as the background one. 



A. Solution for <5 7 

Equation (|3.13|) is easily solved. In this equation, there is no singularity even at the 
background naked singularity. Before the appearance of the background naked singularity, 
i.e., t < 0, we impose a regularity condition to guarantee the locally Minkowskian nature at 
r = 0; 8 1 is C 2 function with <5 7 '| r= o = 0. We also impose the same condition at r = after 
the background naked singularity formation t > 0. Then we obtain general solutions for <5 7 

5 1 = A 7 (i-r) + A 7 (t + r), (4.1) 

where A 7 is an arbitrary C 2 function. Here we should note that after the appearance of the 
naked singularity t > 0, the regularity condition 5 7 '\ r=0 = does no longer guarantee the 
locally Minkowskian nature since the background spacetime itself does not have the locally 
Minkowskian nature at the naked singularity. 



B. Solution for V 

The solution for the Euler equation ()3.17|) is given by 

V = C v (v), (4.2) 

where Cy is an arbitrary function which is everywhere finite. Therefore even at the back- 
ground naked singularity, the velocity perturbation V does not vanish. This fact means that 
the type of the singularity formed at the symmetric center r = is changed from that of 
the background Morgan spacetime. Up to the first order of e, the Ricci scalar is obtained as 

y ,.- 8 , C T V 16 ' Ce ; fl1 ' (4-3) 

We can easily find in the above equation that the Ricci scalar diverges at the symmetric 
center r = if and V do not vanish there. Thus the background naked p.p. curvature 
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singularity becomes the s.p. curvature singularity which is defined as the spacetime singu- 
larity accompanied by the divergence of a scalar polynomial constructed from the Riemann 
tensor 128J. 



C. Solution for 6r 

Using the solution (J4.2j) . we can get the solution for 5r. For notational simplicity, we 
introduce a function 5* defined by 

S(v) = 87rGe tB D B C v . (4.4) 

Note that two times S(v) is the source term in the equation (|3.14j) for Sr. 

As mentioned, we search for the solution which can be expressed by Maclaurin series in 
the neighborhood of the background naked singularity as 

6 R = 5%\t) + 5%\ t )r + 4 2) (t)r 2 + 8$\ty... (4.5) 

The source function S is also assumed to be written by Maclaurin series as 

S(v) = S {0 \t) + S w (t)r + SW(t)r 2 + .... (4.6) 

Substituting the above expressions into Eq. (j3.14|) . we obtain the following equations 

- 25$ = 2S®, 4 0) - 64 2) = 2S^\ 4 1} - 124 3) = 2£« ... (4.7) 

Thus the boundary condition at r = is given by 

6 R '\ r=0 = -S(t). (4.8) 

Since the expansion coefficients S^(t) of the source function are finite, the coefficients 8r 
are also finite even in the neighborhood of the naked singularity. 

The solution for 5r with the boundary condition (|4.8jl is easily obtained by constructing 
a Green function. Since t5r should vanish at r = 0, the Green function Gr for t5r should 
also vanish there. Therefore Gr is given by 

1 f°° dk 

GR(t,r;r,x) = — — sin(fcr) sin(fcc) sm{k(t — r)}. (4.9) 

7T J K 
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Thus the solution for 6r is given by 

2 ft foo 

5 R (t,r) = — / dr i dxG R (t,r;T,x)S(T + x) 
ixr Jti Jo 

+ l{A R (t + r)-A R (t-r)}, (4.10) 
where t\ is some constant, A R is an arbitrary C 2 function. 



D. Solution for ip 

Before the appearance of the background naked singularity t < 0, we impose the regularity 
condition on ip on the symmetric axis r = 0; ip is C 2 function satisfying ip'\r=o — 0. However, 
at the naked singularity < t < v w , this condition is not satisfied because the right hand 
side of Eq. (j3.15|) diverges at r = if S\ r= o ^ 0. Even in such a situation, we can find 
solutions everywhere finite. 

Here we assume that ip can be written in the form of the Maclaurin series around r = 0, 

ijj = + ^\t)r + ^ 2 \t)r 2 + .... (4.11) 

Substituting this expression and Eq. (j4.6J) into Eq. ()3.15|) . we find 

_ ^(D = s<®, ^ (0) - 4^ (2) = - 9^ (3) = S^ 2 \... (4.12) 

Thus the boundary condition at the background naked singularity is given by 

xP'\ r=0 = -S(t). (4.13) 

The above condition guarantees that the solution for ip is everywhere finite. To get the 
numerical values of ip, we need numerical calculations. We will do it in Sec.V. 

Here we should note that the background conical singularity at r = and t > v w remains 
conical one up to the first order. This is because the perturbation variables 5 R , 5 7 and i/j 
behave as those of the regular vacuum spacetime at the background conical singularity by 
the present boundary condition. 



E. Solution for 5^ 



It is easy to integrate Eq. (j3.16|) . Using Eqs. (j4.2|) and ()4.4|) . the solution is obtained as 

"( 2S -c v d 4- s 

2 \ dv 



5 D = i,-5 1 + -\2S- Cy—j^- ) + C D {v), (4.14) 
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where Co is an arbitrary function which is everywhere finite. 



V. GRAVITATIONAL RADIATION FROM NAKED SINGULARITY 

In this section, we study the energy carried by gravitational radiation emitted from the 
naked singularity. Since the energy flux in the vacuum region is given by Eq. (j3.20j) . we 
numerically calculate ip in the causal future of the background naked singularity. In order 
to do so, we adopt the retarded and advanced coordinates u and v instead of t and r. Then 
Eq. (l3.15|) is rewritten in the form 

^ = D u (u,v), (5.1) 
dD v (u,v) dD u (u,v) 1 

7; = i = 7T( T {D v (u, v) - D u (u, v) + S(v)}, 5.2 

ou ov 2{v — u) 

We consider the one-parameter family of the background null dust solutions. We adopt 
the thickness v w of the shell as the parameter and set all the members in this family having 
the same total C-energy E B . Since the total C-energy E B in the background spacetime is 
given by 

E B = \(l- e~ 2 ^) , (5.3) 



where 7^ is the asymptotic value of 7 B for v — > 00. Thus independence of C -energy E-q on 
the thickness v w means that 700 does not depend on v w . We can see in the above equation 
that E-q is less than 1/8. If E-q is equal to 1/8, 7^ is infinite, and then the space is closing 
up in the r-direction. 

For notational convenience, we introduce the dimensionless null coordinates defined by 

x = — and y = —. (5.4) 

v w ' v w 

Then the asymptotic value 700 is written as 

7oo = / dv—— = / dyv w ——. (5.5) 
Jo dv Jo dv 

Since 700 does not depend on the thickness v w of the dust shell, d^/dv should be written 
in the following form, 

^ = -F 7 fe). (5.6) 

dv f w 
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We set the velocity perturbation V being 

V = C v (v) = F v (y), 



(5.7) 



so that the maximal value of V does not depend on the thickness v w . Then we find 

S(v) = = -F,F V =: -F s (y). (5.8) 

We introduce new dimensionless variables defined by 

D x (x,y) = v w D u (u,v) and D y (x, y) = v w D v (u, v). (5.9) 
Then Eqs. ()5.1|) and ()5.2|) take the following forms, 

9 -t = d ir = w^) (D "- D " +Fsh (5 - n) 

Since Fg is the function of y, the thickness v w does disappear in the above equation. This 
means that ip, D x and D y are the functions of x and y but do not depend on the thickness 
v w as long as we use the variables x and y as independent variables. 

A. C-energy 

We calculate C-energy carried by the gravitational radiation from the naked s.p. cur- 
vature singularity formed at r = 0. The background singularity in t > v w is a conical 
singularity at which the source function S vanishes. This means that ip is not generated 
there. Thus we focus on the gravitational radiation generated at the intermediate naked 
singularity < t < v w . For this purpose, it is sufficient to investigate C-energy flux through 
the null hypersurface of v = v w from u = to u = v w (see Fig.l). Using Eq. (j3.20j) . C-energy 
A£ns transfered from the intermediate background p.p. curvature singularity through this 
null hypersurface is given as 

A£ NS = 77^ / du(v w - u)D u 2 (u,v w ) = — / dx(l - x)D x 2 (x, 1). (5.12) 
4G Jo 4G Jo 

From the above equation, we find that ASns does not depend on the thickness v w of the 
dust shell. This fact is different from C-energy carried by gravitational waves generated in 
the causal past of the Cauchy horizon 3] (see also Appendix A). 
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From Eq. ()5.5jh we see that is proportional to F 1 . Since the mean value V m of the 
velocity perturbation V is given by 

1 r v w r 1 
V m = — dvV{v) = / dyF v (y), (5.13) 
v w Jo Jo 

V m is proportional to Fy- This means that D x is proportional to 700 V m since the source 
term F$ is the product of _F 7 and Fy. Rewriting 7^ by the background total C-energy as 

7oo = -~ln(l-8£ B ), (5.14) 

we obtain 

AE NS = Const x {Knln(l - 8E B )} 2 . (5.15) 

Since the absolute value of the source term |i*s| should be much smaller than unity so that 
the high-speed approximation is valid, |Vmln(l — 8Eb)\ should be much less than unity. 
However at present, we can not say whether ASns is large compared to the total C-energy 
E-q or not. We need focus on a specific model of our interest and then calculate 'Const' in 
Eq.dHED. 

B. Example 

Here we consider the following example, 

S (v) = - — V m Hl-8E B )y 2 (l-y) 2 Q(y)Q(l-y), (5.16) 
V(v) = 6V m y(l - y)Q(y)Q(l - y), (5.17) 

where G is Heaviside's step function. Note that the above S and V have been chosen so 
that Eb and V m are the total C-energy of the background spacetime and the mean value of 
V, respectively. We solve Eqs. (|5.1U|) and (J5.ll)) numerically and then estimate Eq. (|5.12J) . 

We set the C-energy of the present system being the same as that of the background. By 
Eq. ()3.18j) . this condition leads to the following condition 01111 = 

5 7 -(r5 R )' = 0. (5.18) 

Eqs.flTIH and (l3~T21) lead to 

d v {5 7 - (r5 R )'} = 8nGe^D B L-^ + 5 D - 2d v (r5 R ) - ^} . (5.19) 
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Therefore in the region of Db > 0, the following equation should hold, 



5 7 - tP + 5 D - 2d v (r5 R ) - - = 0. (5.20) 



ip = <5 7 = Sr = and 5d = — on u = u\. (5-21) 



Here we adopt the following initial data so that Eq. (|5.18|) is satisfied, 

V 
2~ 

It is not so difficult to numerically solve the differential equations (J5.1U)) and (J5.ll)) . We 
set the initial condition if> = and thus D y = at the Cauchy horizon x = 0. This initial 
condition guarantees that there is no ingoing flux across the Cauchy horizon and hence we 
can see the gravitational radiation generated just at the naked singularity. The result of the 
three-digit accuracy is obtained as 

AE NS = 0.0643 x {V m ln(l - 8E B )} 2 . (5.22) 

Here for notational convenience, we introduce a new parameter defined by 

£ = -V m ln(l-8£ B ). (5.23) 

Since the source function S is proportional to this parameter e, e should be much less than 
unity so that the high speed approximation is applicable. Eb is expressed by e and V m as 
E B = (l — e~ £//Vm ) /8 and thus we obtain 

A^ns „ . e 2 



- 0.514 x — . (5.24) 

E B i- e -e/v m ^ > 

We can easily see that the above quantity is monotonically increasing function of both e and 

V m for e > and V m > 0. In the parameter space (e, V m ), e = V m = 1 is the marginal point 

at which the high speed approximation must break down. The upper bound on AE^s/Eb 

for the high-speed collapse is realized at this point, 

max ( —^-) = 0.813. (5.25) 



Eb 

The above value is rather large. The large amount of the gravitational radiation seems to 

be emitted if the collapsing speed 1 — V m is not so large and the background C-energy Eb 

is not so small. On the other hand, in the case of very small e, the emitted energy is 

A V 

— — ~ 0.514 x eV m . (5.26) 
Eb 

Thus the emitted energy will be at most a few percents of the total energy of this system 
in the situation that the present approximation scheme describes the system with sufficient 
accuracy. 
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VI. SUMMARY AND DISCUSSION 



In this paper, we studied the gravitational radiation from the naked singularity formed 
by high-speed gravitational collapse of a cylindrical thick shell composed of dust. In this 
issue, it is very important what boundary condition for metric and matter variables should 
be imposed at the background naked singularity. Here we have adopted the boundary 
condition that all the perturbation variables are expressed in the form of the Maclaurin 
series at the background naked singularity. Then we have found a boundary condition that 
guarantees all the metric and matter variables are finite everywhere at least up to the first 
order approximation. At first glance, it is expected that the spacetime of the high-speed dust 
collapse is not so different from the background Morgan spacetime even in the neighborhood 
of the background naked singularity. The boundary condition adopted in the present paper 
realize what we expected. Further this implies that the high-speed approximation scheme 
is everywhere valid. 

We investigated C-energy sent out by gravitational emission from the naked singularity. 
In the situation that the present approximation scheme is valid, every perturbation variable 
is small and as a result, the emitted energy of gravitational radiation is small. We assumed 
a simple model and found that the upper bound of the emitted energy by the high-speed 
collapse is almost the same as the total C-energy of the system. However, in the case that 
the high-speed approximation describes the system with sufficient accuracy, the emitted 
energy might be at most a few percents of the total C-energy of the system. Thus the 
emission of gravitational waves from cylindrical naked singularity can be gentle process. This 
result seems to imply that the conjecture by Nakamura, Shibata and KN is not necessarily 
valid. However the dust might be unrealistic assumption and hence further investigation 
is necessary. On the other hand, the rather large value of the upper bound on the energy 
of gravitational radiation implies that there remains possibility of the large gravitational 
radiation by the gravitational collapse of a cylindrical thick dust shell with not so large 
collapsing velocity. This is also a future problem. 
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APPENDIX A 



Here we consider the energy sent out by the gravitational radiation generated in the 
causal past of the Cauchy horizon u = 0. The emitted energy AEr from u = «;(< 0) to 
u = is given by 

AE R (v) = ^J u du{v - u)D u 2 (u, v ) = ^j x dx{y - x)D x 2 (x, y), (Al) 

where Xi = U{/v w and y = v/v w is larger than unity so that the estimation of AE R is 
done outside of the shell. In contrast to the case of the gravitational waves from the naked 
singularity in Sec.V, the emitted energy AE R depends on the thickness t> w of the shell. In 
the thin shell limit, v w — > 0, the region of integration in Eq. (jAl|) becomes (— oo,0] since 
Xi — > — oo. Thus if D x (x,y) does not vanish sufficiently rapidly in the limit of x — ► — oo, 
this integral diverges. 

Let us consider the initial data of which ip = at u = it; (thus D v = 0), or equivalently 
ip = D y = at x = Xi. This initial condition guarantees that there is no ingoing energy flux 
across x = x«. Here we assume x ; — 1 and then consider the case of \x\ ^> y. In this case, 
Eq. 1)5.11) ) becomes as 

^^L (Dy . Di + Fs) . (A2) 

Substituting the expression 

D x = A(x,y)e£ (A3) 

into the above equation, we obtain 

f~4<^>-~4(f+4 (A4) 
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Integrating the above equation with respect to y, we obtain 

+ J*F s (z)dzy (A5) 

where the integration constant has been set so that A vanishes in the vacuum region inside 
the shell y < at x = X;. Therefore we find for x ~ x- x and y > 1, 

D x ~ -— , (A6) 

where 

M = \( l F s {z)dz (A7) 

is a positive constant. 

The above result means that if we set x\ being — oo, D x (x,y) — > —M/x for x — > — oo. 
Hence introducing some constant A <C — 1, we see that in the large Xi limit, 



1 r° 

AEr = 4G J dx (y ~ x ) D * ( x >y) 

= {jx + L ) dX< " y -x ) jD ^ 2 ^' ?/ ) 



1 /" A , M 2 
AG 



r A M 

I dx h F(y; A) — > oo for x\ — > — oo, (A8) 

J x\ X 



where F(y; A) is some function of y and the parameter A. This results means that in the 
thin-shell limit with the total C-energy fixed, the emitted energy diverges. This is consistent 
with the result in Ref . [lsj] . 



[1] R. Penrose, Riv. Nuovo Cim. 1 (1969), 252. 

[2] S.W. Hawking, Commun. Math. Phys. 25 (1972), 152. 

[3] P. Joshi, Global Aspects in Gravitation and Cosmology, (Oxford University Press, 1993). 

[4] T. Nakamura, K. Maeda, S. Miyama and M. Sasaki, Proceedings of the 2nd Marcel Grossmann 
Meeting on General Relativity edited by R. RufRni (Amsterdam: North-Holland 1982), p675. 

[5] T. Nakamura and H. Sato, Prog. Theor. Phys. 67 (1982), 1396. 

[6] S.L. Shapiro, and S.A. Teukolsky, Phys. Rev. Lett. 66 (1991), 994. 

[7] T. Nakamura, M. Shibata and K. Nakao Prog, of Theor. Phys. 89 (1993), 821. 



19 



[8] C. Cutler and K.S. Thorne, in Proceedings of the 16th International Conference on General 
Relativity and Gravitation, edited by N.T. Bishop and S.D. Maharaj (World Scientific, 2002), 
p.72. 

[9] H. Tagoshi et al, Phys. Rev D63 (2001) 062001; 

M. Ando, et al, Phys. Rev. Lett. 87 (2001), 3950. 
[10] K.S. Thorne, in Magic Without Magic; John Archibald Wheeler, edited by J. Klauder (Frieman, 

San Francisco, 1972), p231. 
[11] S. A. Hayward, Class. Quantum Grav. 17 (2000), 1749. 
[12] T. Piran, Phys. Rev. Lett. 41 (1978), 1085. 
[13] F. Echeverria, Phys. Rev. D47 (1993), 2271. 
[14] T. Chiba, Prog. Theor. Phys. 95 (1996), 321. 

[15] K. Nakao and Y. Morisawa, Class. Quantum Grav. 21 (2004), 2101. 

[16] K. Nakao and Y. Morisawa, |arXiv:gr-qc"/ 0409042 

[17] T. Harada and K. Nakao, Phys. Rev. D70 (2004), 041501. 

[18] R. M. Wald, General Relativity (The University of Chicago Press, 1984). 

[19] R.W. Wald, J. Math. Phys. 21 (1980), 2802. 

[20] G.T. Horowitz and D. Marlof, Phys. Rev. D52 (1995), 5670. 

[21] A. Ishibashi and A. Hosoya, Phys.Rev. D60 (1999) 104028. 

[22] A. Ishibashi and R.W. Wald, Class. Quant. Grav. 20 (2003) 3815. 

[23] G.W. Gibbons, S.A. Hartnoll and A. Ishibashi, |h~ep-th/ 0409307] 

[24] K.S. Thorne, Phys. Rev. 138 (1965), B251. 

[25] T. A. Morgan, GRG 4 (1973), 273. 

[26] P. S. Letelier and A. Wang, Phys. Rev. D49 (1994), 5105. 
[27] B.C. Nolan, Phys. Rev. D65 (2002), 104006. 

[28] SW. Hawking and G.F.R. Ellis, The large scale structure of space-time (Cambridge University 
Press, 1973), p. 260. 



20 



